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Abstract. Kirillov and Roshotikhin conjectured what is now known as the 
fermionic formula for the decomposition of tensor products of certain finite 
dimensional modules over quantum afBne algebras. This formula can also be 
extended to the case of q-deformations of tensor product multiplicities as re- 
cently conjectured by Hatayama et al.. In its original formulation it is difficult 
to compute the fermionic formula efficiently. Kleber found an algorithm for 
the simply-laced algebras which overcomes this problem. We present a method 
which reduces all other cases to the simply-laced case using embeddings of 
affine algebras. This is the fermionic analogue of the virtual crystal construc- 
tion by the authors, which is the realization of crystal graphs for arbitrary 
quantum affine algebras in terms of those of simply-laccd type. 



1. Introduction 

In 1987 Kirillov and Reshetikhin ||2^ conjectured a formula, now known as the 
fermionic formula, for the decomposition of tensor products of certain finite dimen- 
sional representations over an untwisted quantum affine algebra Uq{g) into its Uq{g) 
components, where g is the simple Lie algebra associated with the affine Kac-Moody 
algebra g. The conjecture is motivated by Bethe Ansatz studies. Recently, conjec- 
tures for fermionic formulas have been extended to q-deformations of tensor product 
multiplicities 0, In type An'' this g-tensor multiplicity formula appeared in psf . 
For a single tensor factor, the fermionic formula gives the g-isotypical components 
of a C/g(g)-module associated with a multiple of a fundamental weight. This conjec- 



ture was proven by Chari g in a number of cases. Recently, Nakajima |27 showed 
in the simply-laced case that the characters of such modules satisfy a certain system 
of algebraic relations (Q-system). Combining the result of [Q, his result completes 
the proof of a "weak" version of the q — I fermionic formula in this case. 

The term fermionic formula was coined by the Stony Brook group | p^ , , who 
interpreted fermionic-type formulas for characters and branching functions of con- 
formal field theory models as partition functions of quasiparticle systems with "frac- 
tional" statistics obeying Pauli's exclusion principle. 

Fermionic formulas are g-polynomials or g-series expressed as certain sums of 
products of q-binomial coefficients 
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where ["+^] ^ {q)rn+p/{q)7n{q)p is the ^-binomial coefRcient with (q)„i n"i(l- 
g*), cc({TOj-"'}) is some function of the summation variables m'f'^ and p^"-* is the 
vacancy number (see (5.2)). The summation variables are subject to constraints 
(O). Those sets {m'f^} satisfying (Q are called admissible configurations. From 
(4.1) alone, it is computationally difficult to find the admissible configurations, 
making the evaluation of the fermionic formula intractable. For simply-laced alge- 
bras g, Kleber |l^, |l^ has given an efficient algorithm to determine the admissible 
configurations {mi'}. This algorithm generates a rooted tree with nodes labelled 
by dominant integral weights such that the tree nodes are in bijection with the 
admissible configurations. For non-simply laced algebras, the algorithm fails: some 
admissible nodes cannot be reached. 

One of our goals in this paper is to modify Kleber's algorithm to work in all 
types. This is accomplished by using the well-known natural embeddings of any 
affine algebra into another of simply-laced type |^ : 

^(1) .(2) ,(2)t r,(2) ^ .(1) 

4(2) r(i) ^ n(i) 

It is not hard to express the fermionic formula of the smaller algebra in terms of the 
larger; we call this the virtual fermionic formula. Our algorithm is an adaptation 
of Kleber's algorithm in the simply-laced affine algebra, which trims the tree so 
as not to generate nodes that cannot contribute to the virtual fermionic formula. 
This algorithm succeeds by using some nodes in the larger weight lattice that do 
not correspond to weights in the embedded weight lattice. 

Fermionic formulas denoted M have crystal counterparts. Crystal bases were 
introduced by Kashiwara Q and are bases of C/g(0)-modules in the hmit g ^ 0. Let 
us denote the one-dimensional configuration sums, which are generating functions 
of highest weight elements in tensor products of finite dimensional crystals with 
energy statistics, by X. It was conjectured in that X = M. 

In light of the embeddings of affine algebras (1.1), one might hope that such em- 
beddings also exist for the quantized algebras. Unfortunately they do not. However 
we assert that such embeddings exist for all finite-dimensional afhne crystals, and 
give a construction for them in terms of crystals of simply-laced type. A virtual 
crystal is such a realization of a crystal inside another of possibly different type. 
Perhaps the first instance of a virtual crystal is Kashiwara's embedding of a crystal 
of highest weight A, into that of highest weight fcA where fc is a positive integer [jl^ . 
Extending Baker's work |^], in p9| we conjectured that finite dimensional crystals 

of type Cn'' , A^2l j ^'^'^ ^i+i can be realized in terms of crystals of type ^2^-1 • We 
proved this for crystals associated with single columns (i.e. fundamental weights). 

In this paper we establish the correctness of the virtual crystal approach for 
crystals associated with single rows (that is, multiples of the first fundamental 
weight) for the two infinite families of embeddings. 

The paper is organized as follows. In Section |2| we review the essentials of crystal 
theory. Virtual crystals are introduced in Section || and the characterization and 
validity of virtual crystals associated with single rows is proven. Sections ^ and 
^ review the fermionic formulas conjectured in ^ and the Kleber algorithm, 
respectively, and describe their virtual counterparts. 
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2. Crystals 

2.1. Affine algebras. We adopt the notation of ||]. Let g be a Kac-Moody Lie 
algebra of affine type A^-*, that is, one of the types Aii\n > 1), Bn\n > 3), 
d'\n > 2), Di'\n > 4), Ei'\n = 6,7,8), F^'K G^i\ > 1), 4f (n > 1), 

^2«-i('^ > 2), > 2), ^ or The Dynkin diagram of fl = Aj[^ 

is depicted in Figure |l| (Table Aff 1-3 in pc|]). Its nodes are labelled by the set 
/ = {0, 1, 2 . . . , n}. Let / = A{0}- 

Every affine algebra q has a simple Lie subalgebra g obtained by removing the 
0-node from the Dynkin diagram. This is summarized in the following table: 



(2.1) 
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Let ai, hi, [i G /) be the simple roots, simple coroots, and fundamental weights 
of Q. Let 5 and c denote the generator of imaginary roots and the canonical central 
element, respectively. Recall that 5 — o.iOa s-nd c = X^ie/ ^^^ii where the Kac 

labels ai are the unique set of relatively prime positive integers giving the linear 
dependency of the columns of the Cartan matrix A (that is, A{aQ, . . . , a„)* ~ 0). 
Explicitly, 

(2.2) 
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Figure 1. Dynkin diagrams for . The enumeration of the 
nodes with / = {0, 1, . . . , rt} is specified under or the right si de o f 
the nodes. In addition, the numbers ti (resp. t^) defined in (2.3) 
are attached above the nodes for r = 1 (resp. r > 1) if and only if 
U + 1 (resp. it ^ 1). 



The dual Kac label is the label for the affine Dynkin diagram obtained by 
"reversing the arrows" of the Dynkin diagram of or cquivalcntly, the coefficients 
giving the linear dependency of the rows of the Cartan matrix A. 

Let P = 0^g^ZAa © Z(5 be the weight lattice of g and P+ = 0^g^Z>oAa. 

Similarly, let P — 0^gjZAa be the weight lattice of g, P = 0^gjZ>oAa, Q — 

^^^j'Zaa the root lattice of g and = ^^^j'ZyoOia with simple roots and 

fundamental weights Ua, A^ for a G I. For X, fi G P write At>/^ifA — //G . 
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For i G / let 

(2.3) = max(— T, On ), = max( — , oo). 

a( a. 

The values U are given in Figure 0. We shall only use i^l and for z G /. For a G / 
we have 

= 1 if r = 1, ta = a'^ \ir>l. 
Let (-j-) be the normalized invariant form on P It satisfies 

(2.4) (a,|a,) = ^Ay 
for i^j^I. In particular 

(2.5) {aa\aa)^% 
if tta is a long root. 

2.2. Crystals. The quantized universal enveloping algebra Uq{Q) associated with 
a symmetrizable Kac-Moody Lie algebra g was introduced independently by Drin- 
feld [Q and Jimbo § in their study of two dimensional solvable lattice models 
in statistical mechanics. The parameter q corresponds to the temperature of the 
underlying model. Kashiwara showed that at zero temperature or g = the 
representations of Uq{Q) have bases, which he coined crystal bases, with a beautiful 
combinatorial structure and favorable properties such as uniqueness and stability 
under tensor products. 

Let g' be the derived subalgebra of g. Denote the corresponding quantized 
universal enveloping algebras of g D g' D g by C/q(g) D C/g(g) 3 C/g(g). 

In 10, ^ it is conjectured that there is a family of finite-dimensional irreducible 
t/^(g)-modules {W^""^ | a G /,« G Z>o} which, unlike most finite-dimensional ?7^(g)- 
modules, have crystal bases _B°''. This family is conjecturally characterized in 
several different ways: 

1. Its characters form the unique solutions of a system of quadratic relations 
(the Q-system) ||. 

2. Every crystal graph of an irreducible integrable finite-dimensional Ug{g)- 
module, is a tensor product of the B""'^. 

3. For A G P let V^(A) be the universal extremal weight module defined in 
Section 3] and B{X) its crystal base, with unique vector u\ G B{X) of weight 
A. Then the affinization of i?"'* (in the sense of ^^) is isomorphic to the 
connected component of ux in B{X), for the weight A = iAq. 

In light of point (2) above, we consider the category of crystal graphs given by 
tensor products of the crystals B'^-^. 

We introduce notation for tensor products of i?"^*. Let 

(2.6) B= (g) (B'^.»)®i^°\ 

(a,j)e7xZ>o 

where only finitely many lI""* are nonzero. In type An"^ this is the tensor prod- 
uct of modules, which, when restricted to A„, are irreducible modules indexed by 
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rectangular partitions. The set of classically restricted paths (or classical highest 
weight vectors) in B of weight A G P is by definition 

■p[B, A) = {6 e B I wt(fe) = A and ab undefined for aU i e 7}. 

Here is given by the crystal graph. For b,b' E B°"^ we have b' ~ ei{b) if there 

is an arrow b' — ^ b in the crystal graph; if no such arrow exists then ei{b) is 

undefined. Similarly, b' — fi{b) if there is an arrow b b' in the crystal graph; 
if no such arrow exists then fi{b) is undefined. If Bi and B2 are crystals, then for 
&i ® 62 e -Bi ® B2 the action of is defined as 



ei{bi ® 62) 



eibi®b2 if ei(&i) > (y£'i(&2), 
61 (g) ei&2 else, 



where ei{b) = max{fc | e^b is defined} and ipi{b) = max{fc | f^b is defined}. This is 
the opposite of the notation used by Kashiwara fill]. 



2.3. Simple crystals. Let W be the Weyl group of g, {si \ i E 1} the simple 
reflections in W. Let B be the crystal graph of an integrable J7g(g)-module. Say 
that & G i? is an extremal vector of weight A G P provided that wt(&) — A and 
there exists a family of elements {bw \ w G W} C B such that 

1. bio — b for w = e. 

2. If {h, , wA) > then e,(6„) = and /i''" ""^^ (&t„) = b.^^^. 

3. If {hi , wX) < then /i(6u.) ~ and e^'^' '"''^^(&u,) — bg.^. 
Following ||l|, say that a C/^(0)-crystal B is simple if 

1. i? is the crystal base of a finite dimensional integrable [/q(g)-module. 

2. There is a weight A G such that B has a unique vector (denoted u{B)) 
of weight A, and the weight of any extremal vector of B is contained in WX 
where W is the Weyl group of g. 

In the definition of simple crystal in condition 1 is not present. However we 
always want to assume both conditions, so it is convenient to include condition 1 
in the definition above. 

Theorem 2.1 (|]). 

1. Simple crystals are connected. 

2. The tensor product of simple crystals is simple. 

For the [/^(g)-crystal B, define e, ip : B ^ P hy 

e(b) = e,(6)A, and ^(6) = ^ (/7,(6)A,. 

iel i€l 

Then the level of B is 



(2.7) 



lev(B) = min{(c, e(6)} \ beB}. 
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2.4. Dual crystals. The notion of a dual crystal is given in |T^, Section 7.4]. Let 
-B be a C/g(g)-crystal. Then there is a t7g(0)-crystal denoted obtained from B 
by reversing arrows. That is, i?^ = {6^ | b S B} with 

wt(5^) = -wt(5) 

(2.8) (^,(5^) = e,(6) 

e.(fo^) = (/.(fo))^ 
/.(6^) = (e.(5))^ 



Proposition 2.2. |l^ There is an isomorphism {B2 ® BiY — B^ (g) B2 given by 



(62® 61)^ H^foV^feV 



22, Section 4] that there 



2.5. One dimensional sums. In this section we recall the structure of a J7g(0)- 
crystal as a graded (0)-crystal. The grading is given by the intrinsic energy 
function D : B Z. For b G B, one may define D{b) as the minimum number of 
times Co occurs in a sequence of operators involving e^, fi for i G I and eg, leading 
from u(B) to b. However we prefer to work with the following concrete definition 
when _B is a tensor product of crystals of the form B^'". This definition essentially 
comes from but it is useful to formulate it as follows [29|. 

Let -Bi,i?2 be simple [/^(b)- crystals. It was shown in 
is a unique isomorphism of C/^(0)-crystals R = Rb2,Bi : B2® Bi Bi® B2, called 
the combinatorial R matrix. In addition there exists a function H : Bi ® B2 ^ I1 
called the local energy function, that is unique up to a global additive constant, 
which is constant on / components and satisfies for all 62 £ B2 and 61 £ Bi with 
i?(&2 61) = 6i b'2 

(2.9) 

if £0(62) > Vo{bi) and eo(&i) > (,00(62) 
H(eo(&2 ® hi)) = H{b2 ®bi)^ {1 if eo(&2) < M^i) and eo(6i) < ^0(^2) 

otherwise. 

We shall normalize the local energy function by the condition H {u{B2)®u{Bi)) = 0. 
It was conjectured in M that 




(2.10) (p{b^) = lev(B''^")Ao for a unique b^ G B 

For a given crystal B^''', denote this element also by u^B^'"). Define the function 
Dgr.s : B''^" ^ Z by 

(2.11) DBr.s (6) = H{b 6") - H{u{B'''') ® b^) 

where H = iJ^r.^^r.., is the local energy function. In all cases in which the Ulj{g)- 
module wi'-* and its crystal base B^''' have been constructed, (2.10) holds and 



(2.11) agrees with the explicit grading on B^''^ specified in a case-by-case manner 
in the appendices of ||^, ||] . 

A graded simple crystal {B, D) is a simple crystal B together with a function 
D : B ^ "L. Let {Bj,Dj) be a graded simple J7^(0)-crystal and Uj = u{Bj), for 
I < j < L. Let B — Bl i?) ■ ■ ■ (i) Bi. Following p8] define the energy function 
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Eb : B ^Ihy 

(2.12) Eb — ^ ifii?i+ii?i+2 • • • 

l<i<j<L 

where Ri is the combinatorial i?-matrix and Hi is the local energy function, where 
the subscript i indicates that the operators act on the z-th and {i + l)-st tensor 
factors from the right. This given, define Db, D'g : B ^ Z by 

L 

D'b = Eb + > ^ DjRiR2 ■ ■ ■ Rj-i 

(2.13) ^tt 

DB{b) ^ D'Bib) - D'bHB)) 

where Dj : Bj ^ Z, acts on the rightmost tensor factor. Then wc say that the 
graded simple crystal {B,Db) is the tensor product of the graded simple crystals 

Theorem 2.3. p9| Graded simple Ulj{g) -crystals form a tensor category. 

Now suppose that for all j, Bj has the form i?'"'" and Dj — Db^ => is the intrinsic 
energy as defined above. Then the function Db is called the intrinsic energy of B. 
Let }}j S Bj be as in ( 2.1C ). Then conjecturaUy there is an element b^ £ B such 

that b^j is the leftmost tensor factor in Rl-i ■ ■ ■ Rj+iRjbK 

Using the Yang-Baxter equation for R and the fact that Rb®b is the identity 
for any B, it follows that 1^ 

(2.14) DB{b) = H{b(g)b^) - H{u{B)(g)b^). 

The one-dimensional sum X{B,X;q) £ ^[q^q^^] is the generating function of 
paths graded by the intrinsic energy 

(2.15) X{B,X-q)= l"""^'^- 

beV{B,\} 

2.6. Crystals of type Bn,Cn, Dn- In this and the next section we will describe 
the classical highest weight crystals B{sAi) and the finite dimensional affine crystals 
J5^'^ for all nonexceptional types as weakly increasing words b in an alphabet X. 
They are also determined by x{b) = (xi)ii£x where Xi is the number of i's in b. 
Whenever an operation yields a negative value for an Xi it will be undefined. 
According to the crystal of B{Ai) has underlying set 

A:" = {1 < 2 < • • • < 71 < o < n < • • • < 2 < T} for 

A:" = {1< 2 < • • • < ?i< n < • • • < 2 < 1} for C„ 

77 — — 

A- = {1< 2 <•••<_<•••< 2 < 1} for L>„. 

n 

The crystal B{sAi) is the set of weakly increasing words of length s in the alphabet 
X such that, in addition, for type Bn there is at most one o, and in type _D„, there 
are either no letters n or no letters n. 

The crystal operators on B{s Ai) are given by 



(2.16) 



,Xi + l,Xt+i ~ 1, . . . ,xi) iixt+i>Xt+i 
, Xi+i + l,Xi - 1, . . . ,xi) if Xi+i < Xi+i 
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with the foUowing exceptions: 



(2.17) 



Type Bn- ej) 
Type C„: e„6 
Type Dn- Cn-ib 

Crib 



, Xji ~h 1, Xq 1, Xtt,, . 



, Xi) if Xq — 

, xi) if Xo = 1. 



, Xn ~t~ 1 , Xf, 



1, 



, Xji^i H~ 1, 1, Xji^ • ■ • 1 -^i) if -^n ^ 

, Xtj, , Xji -\- 1 , X-n, — ! 1,... if^Jyj, — 

, Xt^ — I H~ I7 3^715 -^n 1, . . . , Xi) if X72, ^ 

, Xfi H- 1, Xfn — 1 1, . . . , Xi) if Xji — 0. 



2.7. AfRne crystals B^'^. We recall the crystals B^'^ from [M| (and |20] for type 



Cn^^). The affine algebra g has simple Lie subalgebra of type given in (2.1). There 
is an isomorphism of classical crystals 



(2.18) 



B(sAi) for types si'^ , D^P , ^^'j.^ 

0S(,s'Ai) for types 41 

Bis'M) for type d^\4f. 



The crystal operators ei for 1 < i < n are given in subsection 2.6. The operator eo 
is given by 

(2.19) 



Type i3W,i?W,41i:eo5- 



Type41i?il: eo5 : 



TypeC«,4f: eo5 - 



(a;i,.i:2 - 1, ■ • ■ ,a:2,a;i + 1) if X2 > X2 

(xi — 1,X2, ■ ■ ■ ,X2 + l, Xi) if X2 < X2 

(Xi - 1,2:2, ... ,X2,Xi) if Xi > 
(xi, X2, . . . , :r2, ^1 + 1) if Xi < Xi 

{xi -2,X2, ■ ■ ■ , X2,xi) if xi > ^1 + 2 

{Xi - l,X2, ■ ■ ■ , X2,Xi + 1) if Xi = Xi + 1 

(a:i,a;2, . . . ,X2,xi + 2) ifxi<xi. 



3. Virtual crystals 



3.1. Embeddings of affine algebras. As given in (1.1), there are natural inclu- 
sions of the afhne Lie algebras. These embeddings do not carry over to the corre- 
sponding quantum algebras. Nevertheless we expect that such embeddings exist for 
crystals. Note that every affine algebra can be embedded into one of type y4(i),D(i) 
and E^^^ which are the untwisted afline algebras whose canonical simple Lie sub- 



l(2)t ^(2) 



1(1) 



are 



algebra is simply- laced. Crystal embeddings C„ %v42„ , ^2n ' ^n+i 
studied in ||2§. 

Consider one of the embeddings given in (^]^) of an affine algebra with Dynkin 
diagram X into one with diagram Y. We consider a graph automorphism cr of y 
that fixes the node. For type ^2n-ii '^i'^) ='2n — i (mod 2n). For type D^^^^ the 
automorphism interchanges the nodes n and n+1 and fixes all other nodes. There 
is an additional automorphism for type d'^\ namely, the cyclic permutation of the 
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nodes 1,2 and 3. For type Eq ' the automorphism exchanges nodes 1 and 5 and 
nodes 2 and 4. 

Let and be the vertex sets of the diagrams X and Y respectively, /a 
the set of orbits of the action of a on 7^, and l : ^ 1'^ bijection which 
preserves edges and sends to 0. 

Example 3.1. If X is one of C^^^ , A^^^ , A^^^^ , D^^l-^ and Y = A^^n-i ' t^^en i,(0) = 0, 

= {i, 2n — i} for I < i < n and i(n) = n. 
liX = BrP or ^2ti-i ^'^d F = -D^+i, then = z for i < n and (.(n) = {n, n+1}. 
If X is ' or f]^^ and Y = E^^\ then i(0) = 0, t(l) = 1, i(2) = 3, t(3) = {2, 4} 
and t(4) = {1,5}. 

If X is or G^^' and F = \ then t(0) = 0, i(l) = 2 and i(2) = {1, 3, 4}. 

To describe the embedding we endow the bijection t with additional data. For 
each i G I-^ we shall define a multiplication factor that depends on the location 
of i with respect to a distinguished arrow (multiple bond) in X. Removing the 
arrow leaves two connected components. The factor 7j is defined as follows: 

1. Suppose X has a unique arrow. 

(a) Suppose the arrow points towards the component of 0. Then 7, = 1 for 
all i € I^. 

(b) Suppose the arrow points away from the component of 0. Then 7j is the 
order of a for i in the component of and is 1 otherwise. 

2. Suppose X has two arrows, that is. F = ^2n-i- Then 7^ = 1 for 1 < i < n — 1. 
For z G {0, n}, 7i = 2 (which is the order of a) if the arrow incident to i points 
away from it and is 1 otherwise. 

Example 3.2. For X = B^^ and Y = D^^^i we have 7^ = 2 if < i < n - 1 and 
7„ = 1. For X = A2r(-i ^ — have 7i = 1 for all i. 

The embedding : — > of weight lattices is defined by 

*(Af)=7. 

je<.(i) 

As a consequence we have 

*(af) = 7i E "J 
jei-ii) 

*(5^) = ao7o^''- 

3.2. Virtual crystals. Suggested by the embeddings X ^ Y of afiine algebras, 

we wish to realize crystals of type X using crystals of type Y. 

Let y be a F-crystal. We define the virtual crystal operators Ci, fi for i G I-^ as 
the composites of F-crystal operators fj, ej given by 

f^= n 

iGi.(i) 

^ = n 

These are designed to simulate X-crystal operators fi,ei for i G . The type Y 
operators on the right hand side, may be performed in any order, since distinct 
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nodes j, j' G are not adjacent in Y and thus their corresponding raising and 
lowering operators commute. 

A virtual crystal is a pair (V, V) such that: 

1. F is a F-crystal. 

2. V dV is closed under e^, fi for i ^ . 

3. There is an X-crystal B and an X-crystal isomorphism ^ : B V such that 
ei, fi correspond to Ci, fi. 

Sometimes by abuse of notation, V will be referred to as a virtual crystal. 
Let b <eV and i ^ . We say that b is i-aligned if 

1. <pj(6) = 'pjiib) for all e L{i), and similarly for s. 

2. <fj (b) G 7iZ for all j G and similarly for e. 
In this case 

(3.1) ^f{b)^-^J{^{b}) ior J ^ i{z),b e B 

li ■' 

and similarly for e. Say that 6 G F is aligned if it is i-aligned for all i E I-^ and a 
subset V C V is aligned if all its elements are. 

Proposition 3.3. p9[ | Aligned virtual crystals form a tensor category. 

Say that (V, V) is simple if V and V are simple crystals. For the rest of the 
definitions we assume that the virtual crystals are simple and aligned. 
Let {V,V) and {V\V') be virtual crystals. 

Definition-Conjecture 3.4. Define the virtual R-matrix R" : V ®V' ^ V ®V 

as the restriction of the type Y R-matrix R : V (iSi V ^ V' V . 

For this definition to make sense it needs to be shown that R{V ^ V) C V' ®V. 
In this case, let 'i' : B = V and ^i' : B' = V be X-crystal isomorphisms. By the 
uniqueness of the i?-matrix it follows that the diagram 



B(^B' — ^ B'®B 



(3.2) 



V(g)V' ^" > V'(E)V 



commutes. 



Definition 3.5. Define the virtual energy function : V (S) V ^ Z hy 

H^ib^b') = —HY{b®b') 
70 

where Hy -V ®V' ^1.. 

If Definition-Conjecture holds, it follows that 
(3.3) Hx{b®b')^ H%-^[b)®-^'{b')) 

where Hx : B ® B' "L is the energy function. 



Similarly, define ■.V^'La.s 



D%b) = -D^{b). 
7o 
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If ( |2.14| ) and Definition-Conjecture B.4 hold then 

(3.4) Dx{b) = D'^i'fib)) for 6 GB. 

where Dx : B —> li is the intrinsic energy of B. 
Finally, let A G P for the algebra X and 

ViV, X) = {beV \ wt(6) = ^{X) and e,b = undefined for i G 7^}. 

Then let 

x^v,x)= 

Let us define the F-crystal 

F"-'^ = (g) B"^-"-' 

jei.{r) 

except for A^2n ^'^^ r = n in which case V"'" — By^ <Xi By'^ ■ 

Definition 3.6. Let V^'^ be the subset of V^'^ generated from using the 

virtual crystal operators and fi for i G I'^ . 

Conjecture 3.7. 

(VI) The pair (V^'" jV^'") is a simple aligned virtual crystal. 
(V2) There is an isomorphism of X- crystals 

such that a and fi correspond to Ci and fi respectively, for all i G . 
(V3) Let X be a classical dominant weight for X , B a tensor product of X -crystals of 
the form B^'^ , and {V, V) the corresponding tensor product of virtual crystals 
{V''',V''''). Then 

(3.5) X''{V,X) = X{B,X). 



In (29) Conjecture 13 is proved for emheddings ci^\ A^^J^ , A^^JJ , D^^l^ ^ An-i 



and tensor factors of the form B^ 



Theorem 3.8. Conjecture \3. '\ holds when X is of nonexceptional affine type and 
B is a tensor product of crystals of the form B^ '^. 



This theorem is proven in subsection |3^ 

3.3. Virtual crystals V^ " for A^^J^_^,bIP --^ d';^1^. 

Proposition 3.9. For X = B^^ and Y = D^^j^, 

V^'" = {6 G V^ " I Xt,Xt G 2Z for i < n, Xn + Xn e 2Z, Xn+i = 5„+i = 0}. 
Moreover Theorem [j'.i^ holds. 

Proof. The explicit form of V^''^ follows from u{V^''') — and the definitions of 
the virtual crystal operators. It is easy to show that for s = 1 the map B^'^ — * V^'^ 
defined by i i— *■ m and i ^-^ ii for 1 < z < n and o i— > nn, is the desired isomorphism 
for s = 1. Similarly, it is straightforward to show that for s arbitrary, the desired 
isomorphism \1/ : B^-^ V^'^ is given by replacing each letter (which is an element 
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of B^-^) of a word in B^^'^ by the corresponding pair of letters as in the case s 
This proves (VI) and (V2). 

For (V3) we need to check that 

(3.6) D'^i'fib))^ Deib) for 6 G B. 



Since Db is defined in ter ms o f R, H and functions Db^ ^, it suffices to verify ( 2.14 ) 
and Definition-Conjecture |3.4 

The element u^{B^''^) is given explicitly by P. By the explicit computation of 
H : B^'" (g) B^'" Z given in § it follows that (|1|) holds 



To check Definition-Conjecture p.4| we consider the explicit expressions for the 
i?-matrices of types BjP and D^^^ given in Q. From this it suffices to show that 
the images of relations in the plactic monoid of type Bn are relations in the plactic 
monoid of type D„+i [|6|. This is straightforward. □ 

Proposition 3.10. For X ^ A^2n-i and Y = D^^li, 

V^'' = {be V^^' I Xn+l = Xn+l = 0}. 

Moreover Theorem \3.^ holds. 



Proof. The proof is similar to that of Proposition 3.9, In particular the bijection 
^i.s _^ yi.s gjygj^ |-,y leaving a word unchanged. □ 

3.4. Virtual crystals V^''' for ci^\ A^?^ ^ ^i+i An-i- We require some 
preliminaries on crystals of type ^2ri_i. 

Consider Y = A'^^^^^ and X one of dn\ A^^^ , A'-^^^ D^^}^. In all these cases 
yi'* = _B^"~^'* (g) By^ . We introduce the alphabets 

(3.7) 3^=={1<2< ••• <2n} = {2n^ < (2n- I)'' < ■•• < 2"^ < I''}. 

y and y'^ are the sets of elements of By^ and (_By^)^ ^ B^~^'^ respectively. 



e By" ^'^ is the column of height 2n — 1 in the alphabet y with 



the letter i missing. For 1 < i < 2n - 1, fi{{2n + 1 - iY) = (2n - and f,{b) 
is undefined otherwise, /o(l^) = (2»t.)^ and fo{b) is undefined otherwise. In this 
notation, consists of the weakly increasing words of length s in the alphabet 

3^^. For = 61062 G V^'^ , let yi be the number of letters i in 62 and the number 
of letters in 61, for 1 < * < 2n. 

The i?-matrix R : Bh^ (g) B^^^'^ B^^^'^ (g) Bh^ is given by 



(3.8) 



The i?-matrix R : By'' ® — > B^~^''^ (g) i?y* is given by iterating the above 
i?- matrix so that all of the elements of are commuted to the left. The element 
1^ g) 1 commutes with all elements of By^ and _By"^^'^. 

To formulate the next propositions we also need an involution * : B ^ B on 
crystals of type ^2^-1 lE^l Section 3.8]. Given a word u, let u* be the word obtained 
by replacing each letter i by 2ri -I- 1 — i, and reversing the resulting word. Clearly if 
M is a column word then so is u*. If 6 = C1C2 . . .Cg G B^'^ where Cj is a column word 
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for all j, then by definition b* — c* . . . e B^'^ , which is a sequence of column 
words. Under this map the crystal operators transform as follows: 

Mb*) = e„_,(6)* 

e^{b*) - fn-^{br 

wt{b*) — wowt{b). 



Proposition 3.11. For X = C, 



- 4(1) 



(3.9) V'-' = {be V^'^ I 6^* = Rib), min(yi,j/n,niin(y„+i,y,^+i) e 2Z} 
Moreover Theorem \3.^ holds. 



Proof. We first prove (3.9). By the definition of V^'"^, it suffices to show that the 
right hand side V' of (3.9) contains it(F^''*), and every element of V' is reachable 
from u{V^''') using the virtual crystal operators e^, fi for i e I'^. 
We first digress on the self-duality condition 



(3.10) 



=i?(6). 



By the proof of |2^, Prop 6.8], in the set V^'^, the condition ( j.lOj ) is preserved 
under eo, e„, for 1 < i < n — 1, and similarly for /. For b e V^'", using (| 
equation ( 3.10| ) is equivalent to 



(3.11) 
(3.12) 



y2n+i-j ^ vt - min{y^,y^^) -f min(y,+i, y,^+i) 

vln+i-t = yt- ^Myt^Vi) + ^Myt+i^Vz+i) 



for 1 < i < 2n, where 2/2n+i = yi and i/^ 



n-l-l 



yr- 



We deduce two consequences of (3.1C). Subtracting (3.11) and (3.12) we obtain 



(3.13) 



yi 



y2n+l-i =yi + y2n+l-i 



for 1 < i < 2n. We also have eo(^) = 2/i + 2/2n ~ i^in(2/2n, 2/2ri)- By ( 3.12 ) with 
i = 2n and (3.13) with i = 1, we have 

(3.14) eo(6) -2yi-min(2/i,yJ^). 

Now wc show that u = u{V^''') G V'. This element satisfies yi — sSi^i and 



y'^ = sSi,2n for 1 < i < 2n. Comparing this with (3.11) and (3.12) it follows that u 

satisfies (3.10). It follows that u ^ V'. 

We next check that V aligned. Let b € V and i ^ I-^ . Since b sati sfies (3.10 ) it 
is i-aligned ifl<j<n — Iby Prop. 6.9]. For 0-alignedness, by ( 3.14 ) we see 



that £o(^) is even since min(?/i, yj^) is. The proof that (po{b) is even is similar. So 
b is 0-aligned. The proof that b is n-aligned, is similar as well. So V is aligned. 

Next it is shown that the set V' is closed under ei and fi for i € . Let 
b S v. Bib is self-dual since b is. Note that the quantity min(?/i,?/^) is unchanged 
for i ^ {0, 1}. We have £1(61) — yX and (^51(62) = yi. Hence by the tensor product 
rule, To\n{yi,yi) remains the same upon applying ei. Let z = 0. Since b S V' , b 
is 0-aligned, so that eo(^) G 2Z. Since eo{eob) = eo{b) — 2 is even, by ( 3.14 ), the 
self-dual element eob has the property that min(j/i, j/J^) G 2Z. Thus Cib satisfies 
that property for all i. The property that min(y„_|_i, j/^^j^) G 2Z is satisfied for Cib 
is similar. Thus e^b G V for aU i € . The proof that fib G V for aU i € is 
again similar. 
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Let b ^ v. It suffices to find a sequence of operators Ci and fi leading from b to u. 
We shall induct on the quantity min(2/i, y^), which is invariant under and fi for 
i G /'''^\{0} by previous arguments. Suppose first that e^/ (6) > for some j' ^ 0. By 
alignedness it follows that we may apply a sequence of operators Ci for i G /^\{0} 
to 6, thereby passing to a classical highest weight vector of V^'^ . The classical 
highest weight vectors of V^'^ are given explicitly by Uk = {2n^Y~^\'^^ (g) 1*, for 
Q < k < s. life satisfies min(yi,j/^) = k. By assumption b = Uk for k even. If k — 
then b = uq = u and we are done. If fc > then fob satisfies min(j/i,?/^) — k ~2, 
which is even. We are done by induction. 



We have shown that (3.9) holds and that V^''' is aligned. 

The bijection "if : S^'^ — > V^ '^ is given as follows. Let b G B^'". In the case s = 1, 
the map B^'^ V^'^ is given by z i— > (2n + 1 — i)^ i and i i-^ iy ® (2n + 1 — i). 
The map 5* : i?^'" ^ V^''^ is given by the composite map 

(3.15) B^^' ^ ^ ^ ^ (^2«~i,i)®. ^ (S^^i)®^ 

It follows from ( |3.8| ) that the image of this map is contained in By*. 
Computing this commutation explicitly and using the notation Xi , Xi to describe b 
for 1 < i < ri, and yi,y'^ for 5* (6), we have 

n 

yi = xi - min(a;i,xi) + s - ^(xi + Xt) 

4=1 

n 

(3.16) Ui = Xi — min(a;i, xi) + s — ^^(^i + ^i) 

i=l 

Hi = Xi — min{xi,Xi) + min(a;i_i, Xi_i) for j > 1 
j/j^ = Xi — min(a;i, Xi) + min(a;i_i, x^-i) for i > 1 



To recover yi and for n + I < i < 2n one may use ( p. 11 ) and ( |3.12 ), plus the 



fact that the total number of letters in either bi or 62, is s. 

The composite map given in (3.15) sends e; to for 1 < i < n |p| . It is 



straightforward to check that eo goes to cq using ( 3.16| ), (3.11), and (3.12). There- 



fore 4" is a morphism of X-crystals. It is clearly injective. The image is V^''^ 
since ^'(u(i?^''*)) — u{V^''') and both B^''^ and V^'" are connected. Therefore 
vj; ; jg j^jj^ isomorphism of X-crystals. This completes the proof of (VI) 

and (V2). (V3) follows by ||, Section 6.6]. □ 

Proposition 3.12. For X A^^}^ and Y = 

(3.17) yi-* = {& G yi'* I b^* = R{b), min(2/„+i, y:^+i) G 2Z}. 
Moreover Theorem \3.^ holds. 

The proof is entirely similar to that of 
Proposition 3.13. For X = D^^l^ and Y = A^i^^.i, 

(3.18) V^'" = {6 G V^ " I fe^* = i?(&)}. 
Moreover Theorem holds. 

Proof. For X = -D^+i most of the proof is similar to that of type Here the 

classical subalgebra of X is of type so the isomorphism 5* : B^^'^ V^''^ is a 
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bit different. It is given by o i-> (n + 1)^ (n + 1), with the other letters mapped 
as in the cl%^ case. The exphcit map is given as in (3.16) except that 



(3.19) 



yi = xi — min(xi, xi) + s — Xq — ^^(a;^ + Xi) 

1=1 

n 

uX = xi - min(xi, xi) + s — Xo - ^(a;* + Xi). 

1=1 



□ 



Proposition 3.14. For X = A''^}^'^ and Y = A^^^ 



(3.20) V^^' = {6 e V^^' I b""* = R{b)M^{vi:VX) e 2Z}. 

Moreover Theorem \3.^ holds. 
The proof is similar. 

4. FeRMIONIC FORMULA 

4.1. Review. This subsection reviews definitions of 0|. For this section we 

assmne that g 7^ ^2ri^! fo'" ^^^^ type we refer the reader to Section 7.6]. Fix 

X £ P and B a tensor product of crystals of the form B^'^. Let L^"'' be the number 

of tensor factors in B that are equal to B"'*. Set aa — eta for all a e / except for 
(2) 

type in which case eta are the simple roots of type i?„. 

Let v = {m'f'^) be a matrix of nonnegative integers for i G Z>o and a £ I. Say 
that 1/ is a (B, A)-configuration if 

(4.1) ^ ^m^^aa^ ^ zLi-^^A^-A 

iGZ>o i£Z>o 

except for type A^,^ . In this case the right hand side should be replaced by t(r.h.s) 
where i is a Z-linear map from the weight lattice of type C„ to the weight lattice 
of type Bn such that 

— c, for 1 < a < 



2Af for a = n. 



Say that a configuration is admissible if 



(4.2) pI"^ > for all a e / and i e Z 



where 

(4.3) pi'^) 



j ii"'' min(?, k) - ^(aajab) min(4i, iafc) m-i''^ j 



feez> 

Write C(i?, A) for the set of admissible [B, A)-configurations. Define 



(4.4) cc(i') = I ^ ^ {aa\ab)TiiiTL{tbj,tak)rri^°'''rrii^ 



[a] (b) 
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The fermionic formula is defined by 
(4.5) 



r (a) 



h m. 



(a). 



Tlie X = Af conjecture of [| § states that 
(4.6) X{B,X;q-')=M{B,X;q). 

The fermionic formula M{B, A) can be interpreted using combinatorial objects 
called rigged configurations. Denote by {v, J) a pair where v = {m[°''^) is a matrix 
and J = (J^"'*^) is a matrix of partitions with a & I and i € Z>o. Then a rigged 
configuration is a pair {v, J) such that v € C{B, A) and the partition J^"'') is 
contained in a m'f \h') x p^\v) rectangle for all a,i. The set of rigged {B,X)- 



configurations for fixed A and B is denoted by RC(i3,A). Then (4^) is equivalent 
to 

M{B, A) = J2 9''^""^'' 

(i^,J)eRC(B,A) 

where cc(i., J) = cc(i/) + | J| and \ J\ = E(a.^) *ak^°''^l- 

4.2. Virtual fermionic formula. We define virtual rigged configurations in anal- 
ogy to virtual crystals. 



Definition 4.1. Let X and Y be as in (f .1), and A and B as in subsection 4.1 for 
type X. Let {V,V) be the virtual F-crystal corresponding to B. Then RC^(i?, A) 
is the set of elements (P, J) G RC(V^, ^(A)) such that: 

1. For all i G Z>o, rnf"-* = fh'"''^ and J^"-* — J^'''' if a and 6 are in the same cr-orbit 
in 



2. For all i G Z>o, a G , and G t(a) C , we have mj''-' = if j ^ and 
the parts of are multiples of 7^. 

Theorem 4.2. There is a bijection KC{B, A) RC^(i?, A) sending [v^ J) ^ (9, J) 

given as follows. For all a Q I , b £ t(a) C / ,, and i G Z>o, 



-(b) (q) 



(4.7) 

(4-8) 

(2) 

except when X ~ A^^l and a — n, in which case 

m,- = m] 



The cocharge changes by 
(4.9) 



cciv, J) — 7o cc{v, J). 



Proof. Let L be to F as L is to B as in subsection 4.1. For a G / , 6 G t(a), and 
i G Z>o, 



for i ^ 7aZ, 
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except when X — A'^^ and a — n, in which case L^"' = 2L\"'^ for aU i. Using (4.3) 
we have, for all b G /,(a) and i G Z>o, 

Ab) (a) 

except when X — A'^^ and i ^ n, in which case p^""* = 2p^^\ Therefore {v, J) ^ 



(9, J) defines a bijection. Using (4.4) we see that (4.9) holds. □ 



5. Algorithms for computing the fermionic formula 

To compute the fermionic formula M{B, A), one must find the set of admissible 
(S, A)-configurations C{B,X). One direct approach would be to test the admis- 



sibility conditions ( [4.2| ) on the set of (B, A)-configurations (4.1) which consist of 
all possible n-tuples of partitions of sizes that depend on A and B. This quickly 
becomes infeasible as B and A grow. 

In |l^ Kleber gives an efficient algorithm to compute the set of admissi- 
ble configurations in the simply-laced types An \ Dn \ E^t'\ and -Eg^^ It 
generates a rooted tree T{B) whose nodes are labelled by elements of . The 
tree T{B) is constructed to have the property that the elements of C{B, A) are in 
bijection with the nodes of T{B) labelled A. If a node x labelled A corresponds to 
a configuration i/, then v can be recovered from the unique path in T{B) from x to 
the root. 

5.1. Kleber's algorithm. We review Kleber's algorithm ||l^, |l^. Let X be the 

Dynkin diagram of an untwisted affine Lie algebra whose canonical simple subalge- 



bra is of simply-laced type. Let B and L be as in subsection 4.1. 

We define a tree T{B) by the following algorithm. Each node x is labelled by an 
element wt(2:) G P called its weight. It has the property that if x is a node and 
y is its child, then wt(x) ^ wt(y) and wt(x) > wt(y). A tree edge {x,y) is labelled 
by the element dxy ~ wt(a;) — wt(y) G Q^\{Q}. 

1. Let To be the tree consisting of a single node of weight and set ^ = 0. 

2. Add 1 to ^. _ 

3. Let T'g be obtained from Tt-i by adding Yl,"a=i Yl,i>i ^i^^ the weight of 
each node. 

4. Let Ti be obtained from as follows. Let x be a node at depth ^ — 1 of 
weight /i. Suppose there is a weight t G P^ such that n t, jj, \> t, and if x 
is not the root, v — 2fi + t d Q where v is the weight of the parent w of x. 
In every such case we attach to a: a child y of weight r. Note that if x is not 
the root, the condition — 2/i -I- r G is equivalent to d^x ^ d^y 

5. UTe^ Tf-i then go to step |. 

6. Otherwise set T{B) = Ti and stop. 

For large £ step || does not change the tree. For such i, step ^ can only be applied 
finitely many times since there are finitely many elements of P dominated by a 
given element of P^ . Hence the algorithm terminates. 

There is a bijection from the nodes of T{B) and the configurations C{B) = 
lj^g-p+ C{B, A) given as follows. Let x be a node at depth p in T{B) of weight A. 

Let A(°),A(1),...,A(P) = A be the weights of the nodes on the path from the root 
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of T{B) to X. Then the configuration i> £ C{B, A) corresponding to x is defined by 



(5.1) m'f^ = (A('-i) - 2A(^) + A('+i) | A,) 



where we make the convention that A = A^^^"'^-' 
numbers are given by 



The vacancy 



(5.2) 



j>i 



-Fr+ 



Suppose we are only interested in finding C{B,X) for a particular A G P . It 
is wasteful to generate the entire tree T{B) and then select the nodes of weight A. 
Because the weight of a node dominates that of any of its children, we can prune 
the tree as follows. In step 0, we only add a node of weight r at depth £ if 



(5.3) 



> A. 



There is another condition under which we can prune. Suppose that in the absence 
of pruning, we would have added a node y of weight r at depth £ in step ^, with 
parent x. Then we do not add y if there is an a such that (r' — A | A^) > and 
{dxy I Aa) = 0. For in this case, the condition in step ^prevents one from reaching 
the weight A as a descendant of t. 

Example 5.1. Let S = 53,2^^2,1^^14^51,1 of type yl^^^ The Kleber algorithm 
produces the tree T{B) given in Figure |2[ The corresponding configurations are 
given in the following diagram, where we represent v as a, sequence of partitions 
with to|"^ rows of length i. The vacancy number p^""^ is placed to the right of 
a row of length i in i/^"^ . 



A 



.(2) 



/3) 



2 Ai + A2 + 2 A3 
2 A2 + 2 A3 
Ai + 3A3 
Ai + A2 + A3 
2A3 

3Ai + A3 
A2 

2Ai 



□ 00 

□ 1 Do 

□ 1 □ 1 Do 

W W 

mo mo Do 

Do Do 

mo mi mo 



□ 1 



m 
mo mo 
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Figure 2. Tree T{B) 



5.2. Virtual Kleber algorithm. Outside of the simply-laced case, Kleber's al- 
gorithm does not directly apply. However we use the embeddings of affine algebras 
into those of simply-laced type, where Kleber's algorithm can be applied. We call 
our method the virtual Kleber algorithm. Let X and Y be as in (LI). Theorem 
4.2 defines a bijection C{B, A) = C^{B, A) where C^{B, A) consists of the configu- 
rations P G C{V, ^'(A)) constrained as in Definition Lf , or equivalently, the 9 such 
that (P, J) G RC^(i?,A) for some J. A naive approach would be to run Kleber's 
algorithm to compute the set C{V,'i>{X)) and then to select the desired subset 
C'"{B,X). A more efficient way is to prune the branches that cannot contain ele- 
ments of C'^'{B, A). This results in a good algorithm to find C'"{B, A) and therefore 
M{B, A) for any affine type. 

More precisely, one only adds the child y to the node x in step |^ at depth £ if: 

1. {wt{y) I aa) — {wt{y) \ ai,) if a and b are in the same cr-orbit of . 

2. If £ — f ^ 7a2, then dy^^ = dxy where w is the parent of x. 



These conditions arc equivalent to those in Definition 4.1 



Let T{B) be the resulting tree. Let 7 = maxa7a. Then there is a bijection 
between C^{B, A), and the set of nodes y of weight A in T{B) that satisfy either of 
the following conditions: 

1. y is at depth £ with £ £ 7Z, or 

2. (dxy I Aa) = for every a such that 1 < 7 = 7^, where x is the parent of y. 
Observe that for £ 7Z, there may be nodes at depth £ in Tg whose weights are 

not in the image of the embedding P-^ =— > , but rather in a superlattice of index 
7. These weights, which cannot appear in the final tree, are necessary as they allow 
the virtual Kleber algorithm to reach all of the desired weights. 

Example 5.2. Let A = C^^\ Y = A'^^\ B = B^'^^B^'^^B^^^ . The virtual Kleber 
algorithm produces the tree T{B) given in Figure |[ The nodes corresponding to 
elements of C^{B,X) are circled. We list the configurations corresponding to the 
circled nodes, ordered by increasing depth and then from left to right. Here we 
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Figure 3. Tree T{B) 



represent u as a sequence of partitions z/^"^ with m-"^ rows of length i. The vacancy 
number p^"^ is placed to the right of a row of length i in z^^"' . 



A 



3Ai + A2 
Ai + 2A2 

3Ai 
Ai + A2 

Ai 



/(I) 



/2) 



□ 

□ 1 

m 1 

n 1 

1 






m 
m 2 




/(3) 





□ 

□ 1 
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